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We investigate the dynamics of cavity cooling of a single ion beyond the Lamb-Dicke regime and 
demonstrate a cooling limit of approximately 50% of the Doppler temperature using a high finesse 
cavity for the first time. We also study the statistical properties of the ion-cavity emission and 
present a model that maps the phonon states to the photon states. With this model, we explain the 
super-Poissonian photon distribution observed in our experiments and propose a method to estimate 
the ion temperature in a real-time manner based on the statistical behavior of the photon emission 
from the cavity. 
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I. INTRODUCTION 

Trapped ions have proven to be a promising system for 
quantum information processing (QIP) applications [T]- 
|3] with all of the experimental requirements having been 
demonstrated. In recent years experimental efforts have 
been made towards establishing entanglement between 
distant traps [IH5]. This remote entanglement serves as 
a building block for quantum networks and provides a 
path to large scale QIP. While ions are good candidates 
for stationary processes due to their long lived internal 
states (e.g. quantum memory), photons are the best car- 
riers of quantum information between physically sepa- 
rated sites [345] . Thus, an ion-photon interface is im- 
portant for the development of large scale QIP. An ideal 
system for such an interface is based on an ion trapped 
within a high finesse cavity [T0l - fT2| . The cavity enhances 
the interaction between the ion and a single photon, and 
enables efficient extraction of emitted photons. Proposed 
applications of trapped ion-cavity systems in QIP include 
quantum repeaters [131 H3J , entanglement of distant ions 
[TSUI?] and quantum logic gates [T5H21j . To date, re- 
markable advancements have been made: single photon 
sources [HI [22] , single ion lasers [23] and ion-photon en- 
tanglement [12] have all been demonstrated with trapped 
ion-cavity setups. 

In addition to QIP applications, a cavity also provides 
other useful functions for a trapped ion system such as 
enhanced photon collection efficiency [53], detection and 
minimization of excess micromotion |25j and a means for 
cooling ions [26j . In particular, a notable property of 
cavity cooling is that it can, in principle, be performed 
without affecting the logical information stored in the 
atomic internal states [2"Tll2"5] . In addition, cavity cooling 
techniques can be effective well below the Doppler limit. 
This is important for heavier ions for which Doppler cool- 
ing typically results in a large mean vibrational state, 
making it difficult for further sideband cooling. Cav- 
ity cooling of neutral atoms has been demonstrated and 
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significant experimental milestones such as ground state 
cooling have been achieved [221 30 . However, for ions, 
cavity experiments have been limited [26j and cooling 
below the Doppler limit is yet to be reported. 

In this paper, we investigate cavity cooling of a sin- 
gle ion beyond the usual Lamb-Dicke (LD) regime in 
which the localization of the atom/ion is much smaller 
than the wavelength of the atomic transition. We demon- 
strate ion cooling to sub-Doppler temperatures which, to 
our knowledge, is the first demonstration of sub-Doppler 
cooling of an ion using a high finesse cavity. Motivated 
by the observation of super-Poissonian behavior in the 
ion-cavity emission, we also study the photon statistics 
of the intra-cavity field. It is shown that the thermal 
state of the ion is reflected in photon emission from the 
cavity. This potentially allows for the determination the 
thermal state of the ion from the statistical properties of 
the cavity emission. 



II. CAVITY COOLING 

Two types of cavity cooling have been discussed in 
the literature, namely cavity Doppler cooling and cav- 
ity sideband cooling [28]. The former is used for cooling 
free particles (atoms/molecules) while the latter is used 
when they are strongly confined in the LD regime 28J. 
Here we operate in an intermediate regime where neither 
description is adequate. 

Cavity sideband cooling of trapped atoms has been ex- 
tensively studied theoretically [2H1 1311 132 ■ The cooling 
models are all derived based on the assumption that the 
atoms are in the LD regime, rj 2 (n) -C 1, where rj is the 
LD parameter. In this case, the single atom coopera- 
tivity C [33], which characterizes the ratio between the 
scattering rate into the cavity versus that into the free- 
space, is unaffected by the change in the temperature of 
the atom. As r} 2 {n) increases, the thermal dependence of 
the scattering rate into the cavity has to be taken into 
consideration. Even though our system has a maximum 
f] 2 (n) of only ~ 0.3, the thermal effect is observable. Here 
we experimentally demonstrate that cavity cooling is well 
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FIG. 1. (Color online) The experimental setup. A sin- 
gle 138 Ba+ ion is trapped and coupled to a high finesse cav- 
ity. The emitted photons from the cavity are collected into 
a fiber-coupled single-photon counting module (SPCM). A 
CCD camera, interchangeable with another free-space SPCM, 
detects the fluorescence of the ion in the direction of the black 
arrow. The state manipulation and the ion cooling are per- 
formed with the D\ — D$, R p ,Ri and R2 laser beams. See 
text for a detailed description of each beam. 




FIG. 2. The transition driven by the cavity probing beam R p 
and the intra-cavity field. The metastable D3/2 state is omit- 
ted in the figure. A c is the detuning of the laser frequency 
from the S1/2 P1/2 transition while 8 C is the relative de- 
tuning between the laser and the cavity resonance. Ql is the 
Rabi rate of beam R p . g is the coupling strength between the 
ion and the cavity field, k is the cavity field decay rate. In 
cavity cooling experiments, the induced transition is resonant 
with the red sideband of the z vibrational state by setting 
5 C = — uj z — —400 kHz. In photon statistics experiments, 5 C 
is set to zero. A c is fixed at —110 MHz for all experiments. 



described by an adaption of the model given in [33] to ac- 
count for the dependence of the scattering rate into the 
cavity on the thermal state of the ion. 



A. Setup 

The experimental setup, illustrated in Fig. [T] consists 
of a high finesse cavity aligned with its optical axis trans- 
verse to a linear Paul trap [341435] . Details of the ion trap 
have been reported elsewhere [37]. Briefly, a 5.3 MHz 
RF potential with an amplitude of 125 V is applied via 
a step-up transformer to two diagonally opposing elec- 
trodes. A small dc voltage applied to the other two 
electrodes ensures a splitting of the transverse trapping 
frequencies. Axial confinement is provided by two axial 
electrodes separated by 2.4 mm and held at 33 V. Us- 
ing this configuration, we achieve trapping frequencies of 
2tt x (1.2, 1.1, 0.40) MHz for a single 138 Ba + ion. 

Doppler cooling is achieved by driving the 6S1/2 — > 
6P1/2 transitions at 493 nm and repumping on the 
5D 3/2 -> 6P 1/2 transitions at 650 nm ^ [33]. The 
493 nm cooling laser (D\) and the 650 nm repumping 
laser (D2) are both red-detuned by ~ 15 MHz for op- 
timum cooling. Both lasers are combined into a single 
optical fiber and sent into the trap along the (z — y) di- 
rection. Additionally, D-y and D2 are linearly polarized 
perpendicular to a magnetic field of approximately 1.5 
Gauss. This configuration avoids unwanted dark states 
in the cooling cycle. 

The dual coated high finesse cavity is approximately 
5 mm long and has a finesse of ~ 8.5 x 10 4 at 493 nm and 
~ 7.5 x 10 4 at 650 nm. The cavity is slightly birefrin- 
gent with polarization modes split by 239(1) kHz. The 



polarization modes happen to be aligned to within a few 
degrees of the y and z axis respectively. The cavity length 
is stabilized via the Pound-Drever-Hall technique [35] to 
an electro-optical modulator (EOM) sideband of a low 
linewidth 650 nm laser [33] . Changing the EOM drive fre- 
quency allows us to tune the cavity resonance relative to 
the fixed frequency of the 650 nm locking laser. The laser 
frequency is referenced to a temperature stabilized zero- 
dur cavity and is approximately 500 GHz detuned from 
the repump transition, thus having no impact on the cav- 
ity dynamics. A probe laser (Rp) at 493 nm drives the 
cavity induced Raman transition as depicted in Fig. [2] 
The laser is sent into the trap along the (y + z) direction 
and linearly polarized along the magnetic field direction. 
Additionally, the laser is red-detuned by « 110 MHz from 
the S1/2 P1/2 transition and referenced to the fixed fre- 
quency of the 650 nm locking laser via a transfer cavity 
to ensure that it has a well defined detuning relative to 
the cavity resonance. 

The cavity QED parameters are (<7oi k j7) — 27r x 
(2.13(3), 0.172(5), 10.05(35)) MHz [3JSD], yielding a sin- 
gle atom cooperativity Cq — ffo/ K 7 = 2.64(14) [33] . 
However, the effective cooperativity is greatly reduced 
from this value due to the following factors: the re- 
duced dipole matrix element for 7r-polarization (0.33), 
cavity birefringence (0.672(13) at the midpoint between 
the modes) and branching ratio for S1/2 *H> P1/2 transi- 
tion (0.756(12)) [39]. The combination of these indepen- 
dent parameters results in an effective cooperativity of 
C = 0.446(21), valid for an ion in the LD regime. 
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FIG. 3. The Raman transition and the state detection scheme: (a) The ion is prepared in Si/ 2 |mj = 1/2, n z — n) using a 
493nm<r + laser beam, D3. Afterwards, the Ri and i? 2 Raman beams drive the A-type Raman transition and transfer the ion 
to D 3 / 2 |mj = 3/2, n z = m). Both Ri and J? 2 are detuned by Ar w —500 GHz from Pi/a- (b) The unsuccessfully transferred 
ion from S1/2 is shelved by the D4 laser beam to D 5 / 2 via P3/2 with 87% transfer efficiency, while in 13% of the events, the ion 
will end up in D 3 / 2 . (c) The Doppler cooling beams, Di and D 2 , are turned on. If the ion is in D 5/ / 2 , it is out of the cooling 
cycle and thus no fluorescence will be detected, (d) The repumping beam, D5, depletes the population in D 5 / 2 and moves the 
ion back to the Doppler cooling cycle. 



B. Ion temperature measurements 

A standard technique for determining the mean vibra- 
tion quanta, (n), for a thermal state is to compare the 
transfer efficiencies when driving red and blue motional 
sidebands of a Raman transition between two internal 
states [41] , In this case the ratio of transfer efficien- 
cies can be directly related to the mean vibrational state. 
However for moderate values of (n), this method becomes 
inaccurate as the ratio saturates to unity. For larger (n) 
experimentalists have used another technique which com- 
pares Rabi flopping on the red and blue motional side- 
bands [12 03] . We have found that a more accurate and 
faster determination can be achieved by comparing the 
transfer efficiencies of several motional sidebands in the 
Raman spectra. As (n) increases higher order motional 
sidebands appear. Thus (n) in the direction of interest 
can be obtained from the fit of the Raman spectrum to 
the temperature-dependent probability of transfer from 
S1/2 to D3/2 levels. In general this transfer probability is 
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where 5 m — S — muj, tr is the Raman pulse length, 5 
is the detuning of the Raman frequency from the carrier 
transition, w/27r is the trap frequency along the direction 
of interest and fijj is the effective Raman Rabi rate [41] 
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In this equation, rjR Z is the LD parameter, L"(x) is the 
generalized Laguerre polynomial and the parameter Q r 
is given by 
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where and i7^j 2 are the Rabi rates of the i?i and R2 
lasers, respectively and w 500 GHz is the detuning of 
the Raman lasers from the ionic resonances, S1/2 O P1/2 
or P 1/2 O D 3/2 . 

The procedure for obtaining the spectrum is illus- 
trated in Fig. [3] Briefly, the ion is first prepared in 
Si/2\ m J = 1/2) by switching on a 493nm a + beam, 
£*3, and the 650 nm repumping beam, D2, for 20 fis. 
Afterwards, the Raman beams, i?i and i?2, are turned 
on for 7 [j,s to transfer the ion from Si/2\mj = 1/2) to 
D 3 /2|wj = 3/2) via a A-type Raman process. The un- 
successfully transferred ion is shelved to D 5 / 2 via P3/2 
by a 455 nm beam, D4, which is switched on for 20 [is 
after the Raman pulse. For state detection, the Doppler 
cooling beams, D\ and D2, are turned on for 800 /is. If 
the ion is in D 5 / 2 , it is out of the cooling cycle and thus 
no fluorescence will be detected. If the ion is in D 3 / 2 , 8 
photons on average will be collected into a single pho- 
ton counting module (SPCM). The state detection ef- 
ficiency is only w 87 % due to the imperfect shelving 
process which relies on the relative branching ratio of 
the spontaneous decays between P 3 / 2 — > D3/2 (13%) and 
P3/2 —> D5/2 (87%). After the detection, a repumping 
beam at 614 nm is turned on for 100 /is to deplete the 
population in D 5 / 2 - To obtain a complete spectrum, the 
same procedure is carried out for a range of Raman fre- 
quencies which cover the relevant vibrational resonances. 
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FIG. 4. Three examples of Raman spectra used for inferring 
the temperature. The dots with error bars are the experimen- 
tal data while the solid lines are the fit of Eq. Q. The zero po- 
sition in the horizontal axis denotes the carrier resonance. Fit- 
ted parameters: (a) (n z ) = 2.5(2) and Q r = 2tt x 113(l)kHz. 
(b) (n z ) = 30(1) and ft r = 2tt x 117(2)kHz. (c) {n z ) = 57(2) 
and Sir = 27r x 118(2)kHz. The reduced \ 2 values of all fits 
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FIG. 5. The average phonon number (n z ) versus cavity cool- 
ing time. The unfilled squares are the experimental data for 
the cavity cooling from the Doppler limit (n z ) w 25 along 
the z axis. The solid line is the numerical fit of Eq. ( |B6[ ) 
to the experimental data (black filled circles) obtained in 
the experiment where the ion phonon number is monitored 
at different cavity cooling duration after the initialization at 
(n z ) w 2. The values of the fitted parameters are C = 0.33(3), 
Qlo = 2tt x 23(1)MHZ, d(n r )/dt = 6.3(1.6)/ms and the ini- 
tial temperature (n z o) = 2.4(2). The reduced \ 2 °f the fit 
is w 1.8. The dashed line is the simulation based on the fit- 
ted parameters from the same experiment without considering 
the recoil heating in the radial axes. The dotted lines is the 
simulation of the cavity cooling from the Doppler limit using 
the same parameters. The error bars on the data points are 
the standard deviations obtained from the temperature fits 
similar to those in Fig. [4] 



Examples of the obtained Raman spectra are shown in 
Fig. [4] 

C. Experiments 

In order to drive the cooling sideband along the z axis, 
the cavity resonance is red-detuned by the axial trapping 
frequency u> z ~ 400 kHz from the ion-cavity Raman res- 
onance. To investigate the cooling dynamics, the ion is 
first cooled to the Doppler temperature using a Doppler 
cooling pulse of 1.5 ms. The ion is then cavity cooled for 
a time r and then the temperature is probed. The re- 
sults for a range of times r are presented in Fig. [5] (empty 
squares). A minimum temperature of (n z ) = 13.6(6), or 
54(3)% of the Doppler limit, is achieved by applying a 
cooling pulse of 2 ms. 

To investigate the cooling dynamics further, another 
set of experiments was performed in which the ion was 
initially prepared with a low (n z ) using Raman sideband 
cooling [33] immediately after the Doppler cooling stage. 
In our system, the sideband cooling involves a two-color 



Raman transition used for the ion temperature measure- 
ment, in a setup similar to [36] . The Raman beams, 
Ri and i?2, were arranged in such a way that the state 
transfer was predominantly sensitive to motion along the 
z axis. This was achieved by having a much larger LD 
parameter along the z axis, r\R Z — 0.15, than that along 
other axes, r]R X « r]R y = 0.01. After 3 ms of sideband 
cooling, the ion is then exposed to cavity cooling and 
the temperature probed as before. The resulting (n z ) 
are shown in Fig. [5] using black filled circles. As demon- 
strated, the time evolution of the ion temperature is a 
combination of an exponential relaxation to a steady- 
state temperature and a continuous, slow linear increase 
over time. 

In order to account for our experimental results, we 
adapted the model given in [32] to include for effects 
beyond the usual Lambe-Dicke (LD) limit. A detailed 
account of our model is given in Appendix [X] and Ap- 
pendix [B] Essentially, we consider the effect of the ion 
temperature on the scattering rate into the cavity and 
derive an expression for an effective single atom cooper- 
ativity as a function of (n) along each dimension. This 
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effective cooperativity is substituted into the rate equa- 
tions quoted from [32] , which are subsequently used to fit 
the experimental data. The linear increase in (n) is due 
to recoil heating in the uncooled dimensions. This re- 
sults in a decreasing cooperativity and thus a diminished 
cooling rate. 

In Fig. [5j the solid line is the numerical fit using the 
rate equations Eq. (B4 B5 B6) for a birefringent cavity. 



The fitted value of the cooperativity extrapolated to the 
LD regime, C — 0.33(3), is only 74(7)% of the estimated 
value, 0.446(21). The discrepancy could be due to the 
additional heating along the ion radial direction such as 
recoil heating during the sideband cooling and the envi- 
ronmental heating caused by electronic noise. The pho- 
ton recoil contributes a heating rate of rs 0.05 phonon 
per sideband cooling cycle. In 3 ms of sideband cooling 
process, there are 120 cooling cycles in total. Thus the 
phonon occupation number can increase by ~ 6 phonons. 
The environmental heating is measured to be ss 1 phonon 
per ms. Hence, the sum of these heating effects should 
raise the radial vibrational occupation number by « 9 in 
3 ms which causes an underestimation of the fitted co- 
operativity by 13%. Including these influencing factors, 
the fitted value of the cooperativity extrapolated to the 
LD regime is C — 0.38(3), within 2a range of the esti- 
mated cooperativity, 0.446(21). The rest of the fitting 
parameters are within their respective estimations based 
on the independent measurement of the laser power and 
the sideband cooling efficiency. Using these fitting pa- 
rameters, we simulate the cooling for an ion initialized at 
the Doppler limit and the result is shown in Fig. [5] as a 
dotted line. The dashed line is the result of a simulation 
ignoring the recoil heating in the radial directions. From 
this we conclude that this heating does not significantly 
limit the cooling. Hence, to get an even lower (n z ), a 
higher single atom cooperativity is required. This is pos- 
sible only by improving the experimental setup, such as 
having a cavity with a better finesse or without birefrin- 
gence. 



III. PHOTON STATISTICS 



In the course of our investigations of cavity cooling it 
was found that, when the probe laser was tuned near to 
the cavity resonance, the cavity emission became signif- 
icantly non-Poissonian, as illustrated in Fig. [6] In this 
section we show that this observation can be explained 
by accounting for the dependence of the vibrational state 
on the cavity emission. A theoretical model of phonon- 
photon coupling under the laser-cavity resonant condi- 
tion is presented. It is shown that the thermal state of 
the ions motion can be inferred from the photon distribu- 
tion and vice- versa and this relationship is demonstrated 
experimentally over a range of thermal states. 
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FIG. 6. A typical ion-cavity emission profile. The solid 
line is the mean of the photon count {n) c while the dashed 
line is the variance a^. al is approximately 1.7 times larger 
than (n) c around the cavity peak, indicating super-Poissonian 
statistics. 



Model 



The photon statistics of the spontaneous emission from 
single atoms/ions driven by cw lasers are known to 
be sub-Poissonian if the detection window (r) is small 
(t -C I/7) and approximately Poissonian if the detection 
window is large (t» 1 /j) [33 HB] . In the case of ion- 
cavity coupling, the emitted photons do not follow the 
same statistical behavior as the coupling is governed by 
the cavity-induced Raman process. In this process, the 
photon state is coupled to the vibrational states of the 
ion which in turn affects the photon statistics. To obtain 
the phonon-photon statistical relationship, the Hamilto- 
nian of the system is derived. Based on the Hamiltonian, 
a master equation for the density operator can be estab- 
lished. The master equation can then be solved to obtain 
a steady state solution for the density operator. Further- 
more, the mean photon number and the variance can 
be computed and the statistical relationship between the 
phonon states and the photon states can be evaluated. 

For a trapped ion beyond LD regime, the spatial exten- 
sion of the ionic wave function representing the center-of- 
mass motion (phonon) is no longer small compared with 
the laser wavelength and results in a coupling between 
the ion motion and the driving field [371 LTH) ■ To get the 
Hamiltonian for a trapped ion interacting simultaneously 
with a intra-cavity field and a laser radiation field, we use 
the result for a Raman-type interaction from [48 which 
describes the coupling between an ion and two Raman 
lasers. The expression is modified for the ion-cavity sys- 
tem by replacing a Raman beam with the intra-cavity 
field. When the cavity is resonant with the laser fre- 
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quency, the interaction Hamiltonian is 

Him = hn BO f(ala x ;fj x ) f(ala y ;f] y ) f(a\a z ;ij z ) 

x (at + c) , (4) 

where the operator- valued function /(ata^-; fjj) is defined 
in its normally ordered form as 



= exp [-rfj/2] 



(5) 



1=0 



In both equations, and are the phonon creation 

and annihilation operators, respectively, along the j axis, 
for j G {x, y,z}, while & and c are the respective cav- 
ity photon creation and annihilation operators. The LD 
parameters, fjj, are given by 



1 1 

Vx = Vx , r)y = 77= Vy , Vz = Vz , 



(6) 



where rjj = ky / 'h/(2mujj), k is the wavenumber and m is 
the mass of the ion. The von Neumann equation for the 
density operator (p) is 



dp 



[H int ,p] +£{p), 



dt h 
where C, is the superoperator, 

C(p) = n(2 c) pc — pc 1 c — c'cp) . 



(7) 



(8) 



Since \n) is an eigenstate of f(a^a;fj) with eigenvalue 
exp [— rj 2 /2] L n (rj 2 ), the steady state solutions of Eq. Q 
are of the form 



|n x )(n. x | ® \n y )(n y \ (g) \n z )(n z 



(9) 



where \rij) are phonon Fock state along j axis and 
|ci!n ;E ,n H ,n i ) are photon coherent states with amplitude 
given by 
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We assume a thermal distribution of phonons and so the 
reduced steady state density operator for the intra-cavity 
field is given by 



Pss ^ ^ Plmn |^/,m,n) (^^,m,n I 
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with 
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The average intra-cavity photon number, (n) c , is then 
given by 



(ri) c = Tr{p ss c f c} 
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and similarly 



= Tr{p ss (etc) 2 } 



^ Plmn(\airnn\ 2 + \ai 'I; 



(14) 



With Eq. (13) and Eq. (14), it is straight forward to 
compute the variance u 2 = (n 2 ) c — (n) 2 . 

The photon statistics can be characterized by the Fano 
factor F [49] which is the ratio of variance a 2 to the mean 
(n). In our case this is given by 



F 



(n) 



= 1 
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where all summations are over the indices /, m and n. 

The photon distribution for each individual n- 
vibrational state is Poissonian as given by the properties 
of coherent states. Thus, the effective photon distribu- 
tion is given by the sum of Poissonians weighted accord- 
ing to the thermal distribution. When the ion is cold and 
the phonon population is distributed only around the 
motional ground state, the photon distribution is near 
Poissonian. As (n) increases, F increases indicating a 
super-Poissonian distribution. 



B. Experiments 

To obtain the Fano factor as a function of (n z ), we 
performed experiments on a single 138 Ba + ion prepared 
at several values of (n z ). The control of (n z ) was achieved 
by varying the parameters of the Doppler cooling lasers. 
For (n z ) smaller than the Doppler limit ((n z ) « 25), 
Raman sideband cooling was used. 

In each experimental cycle, a Doppler cooling pulse was 
turned on for 1 ms. This pulse was followed by 3 ms of 
Raman sideband cooling for the measurement with (n 2 ) 
smaller than the Doppler limit. The ion-cavity emission 
was then observed while probing with beam R p . To avoid 
heating of the ion during probing, this beam was switched 
on for only 100 /is. During each probing event, the emit- 
ted photons were counted with a single photon counting 
module (SPCM) coupled to the output of the cavity via 
a single mode fiber. From the transmission of the cavity 
at 493 nm (24%), the fibre coupling efficiency (70%), and 
the quantum efficiency of the SPCM at 493 nm (45%) 
we estimate an overall detection efficiency of intra-cavity 
photons of approximately 7.5%. 
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FIG. 7. Fano factor F as a function of the averaged phonon 
number (n z ). The solid line is the fit of the theoretical model 
to the experimental data while the dashed line simulates the 
ideal situation where the temperatures of the ion are equiv- 
alent in all dimensions. The F error bars are the standard 
deviation of the Fano factor over 300 sets of photon statistics 
while the (n z ) error bars are the uncertainties from the tem- 
perature measurements. Raman sideband cooling was used 
for the first data point such that the (n z ) is brought down 
to ss 2. Note that the sideband cooling is only along z axis, 
hence the ion temperature in other axes should be higher than 
the Doppler limit. This is a one-parameter fit with the fitted 
value of Q R0 = 2tt x 0.14(1) MHZ. The reduced \ 2 of the fit 
is » 0.8. 



Each statistical photon distribution was obtained from 
1000 probing events, yielding a single value for the Fano 
factor. 300 sets of these distributions were collected in 
order to acquire sufficient samples of mean and variance 
values to compute the dispersion of the Fano factors for 
each (n z ). 



C. Results 

Fitting the experimental data to the theoretical model 
requires the information on phonon occupation numbers 
of all vibrational modes. Due to the current configura- 
tion (see Sec. II C I, sideband cooling and temperature 



measurements can be performed only along the z axis. 
However, the ion is assumed to be in thermal equilibrium 
after Doppler cooling so that the temperature along the 
other axes is also at the Doppler limit. We therefore use 
the relationship 



(16) 



which reduces the independent variables to only (n z ) 
and fljiQ. In the measurement where sideband cooling is 
used, the temperature along the non-cavity cooled axes 



is higher than the Doppler limit. However, the additional 
heating can be properly accounted using the approach in 
SecHTUl 

Both the experimental data and the fit are shown in 
Fig. [7] Eq. (151 is fitted to the experimental data where 
the only free parameter is the ion-cavity coupling rate 
fi.R0. The fitted value is 2n x 0.14(1) MHz, which corre- 
sponds to a laser Rabi rate of f2^ = 27r x 20(1) MHz and 
agrees with an estimation based on the measured power 
of the R p laser. Minimizing the x 2 of the fit produces a 
reduced \ 2 of 0.8, indicating that the experimental out- 
come is in good agreement with the model prediction. 

The experimental results confirm that the Fano fac- 
tor is related to the ion temperature as described by the 
model. However, as seen in Fig. [7j the span of the error 
bars reduces the precision in deducing the ion tempera- 
ture for a given Fano number. The Fano error bars in the 
figure are mainly due to the intensity fluctuation of ±5% 
at the cavity output. This is caused by the frequency 
drift in the Doppler cooling laser and the drift of the ion 
position from the cavity anti-node throughout the exper- 
imental period, which slightly affect the ion temperature 
as well as the ion-cavity coupling efficiency. Although 
the error bars are only a few percent compared to the 
signal amplitude, the error bars between the adjacent 
data points are still overlapping due to the small gradi- 
ent curve in the figure. A better precision is achievable 
by improving the laser reference stability which results in 
a smaller error bar or having a larger ion-cavity coupling 
rate, f2.Ro> which results in a steeper curve. As increasing 
the laser intensity £Ilo would induce undesired ion heat- 
ing, a larger FIro can only be achieved by improving the 
ion-cavity coupling strength (go). However, this is not 
readily achievable with our current setup. 

The experiment demonstrates the possibility of us- 
ing the photon statistical distribution as a temperature 
probe. The advantage of implementing this technique 
is that the ion temperature can potentially be measured 
within a relatively short period of time. Applying the 
typical measurement methods as discussed in Sec. |IIB| 
or in |42j usually require long experimental times to ac- 
quire sufficient data. However, the method proposed here 
requires only one statistical distribution to estimate the 
Fano factor for a particular photon distribution. Thus, 
the ion temperature can be estimated in almost a real- 
time manner: less than a second in our experiment to ac- 
quire a statistical distribution. This is particularly useful 
for setup and optimization of processes such as Doppler 
cooling in which rapid monitoring of the ion temperature 
is definitely helpful. 



IV. CONCLUSIONS 

We have demonstrated cavity cooling of an ion to sub- 
Dopplcr temperatures for the first time. A theoretical 
model which describes the cooling dynamics and the lim- 
iting factors is presented and validated by the experimen- 



tal results. Useful information, such as the effect of re- 
coil heating on the cooling performance, are provided and 
may help future efforts in achieving ground state cavity 
cooling in an ion-trap. 

In the second part of the paper, we have shown the con- 
nection between the statistical distribution of the cavity 
photon emission and the ion temperature. This poten- 
tially allows for the use of the photon statistics as a fast 
temperature probe for a trapped ion. 



Approximating uj x 



u> r , Qr is then evaluated as 
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Appendix A: Thermal effect on ion-cavity coupling 

To estimate the thermal effect on ion-cavity coupling, 
we consider the setup depicted in Fig. [T] in which a single 
ion is located at r = xx+yy+zz. The ion interacts with 
a far detuned laser radiation field R p and an orthogonal 
cavity field. The ion-cavity Raman rate fi^ due to the 
cavity induced Raman transition is 



and 



x exp 



ik : 



V2 



exp 



y 2 + z 2 

2 



(Al) 



where VL l (?) is the laser Rabi rate with a maximum value 
of f2z,o at the center of the mode waist, g(r) is the ion- 
cavity coupling rate with a maximum value of go at the 
cavity anti-node, A is the detuning of the laser frequency 
from the ionic resonance, k is the wavenumber, ip c is the 
phase of the intra-cavity standing wave and wl is the 
beam waist of the laser. Due to the position spread of 
the ion in the trap, this Raman rate has to be averaged 
over the Gaussian localization of the ion wavepacket 
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(A2) 



where 
P(r) 



Here <jj is defined as the wave function spread of the ion 
along the j axis 



(27T) 


3/2 o x a y 












x exp 


' x 2 ' 

. 2a 2 _ 


exp 


[ y 2 

[ 2< 


exp 


' z 2 ' 
. 2a* _ 


(A3) 



j' 2 

3 



((TjW L ) 2 



Iff 



2a 2 



(A7) 



If we consider that the ion is positioned at the cavity anti- 
node and the laser beam waist is much larger than the 
wave function spread, <j) c = tt/2 and a'^ rs Uj. Eq. (A5| 
can thus be simplified to 



£1r — Qro exp 



3(jfcer r ) 2 (ka z ) 2 



(A8) 



To obtain a Raman rate with explicit dependency on tem- 
perature, aj can be written in terms of the LD parameter 
rjj and the average phonon occupation number (rij) 



ka-j 



lh(( nj ) + 1/2) 



(A9) 
(A10) 



with 



fc v /fi/(2r 



Substituting Eq. (A10) into 



Eq. (|A8| , a Raman rate with thermal effect incorporated 

3rj 2 ((n r 



is obtained 

Mr = 0_ro exp 



1/2) 



x exp 



2 

1/2) 



(All) 



The single atom cooperativity determines the scattering 
rate into the cavity. Since this scattering rate is pro- 
portional to f2|j, we use the following expression for the 
thermally averaged cooperativity 

C «n r ), (n»» = Cexp [ - 3^ 2 ((n r ) + 1/2) ] 

x exp[-ry 2 (K) + l/2)] , (A12) 

which converges to the single atom cooperativity C for 
an ion in the LD regime r] 2 (nj} <Si 1. For a trapped ion 
beyond the LD regime, we use Eq. ( A12 1 for the effective 
cooperativity at a given temperature. 



Appendix B: Recoil heating vs cavity cooling 



k R T 



(A4) To model our results, we first note that the motional 

sidebands are well resolved by the cavity, thus cooling can 



9 



be selectively performed along a single direction. Other 
directions are then only affected by recoil heating. For 
cooling along a particular direction, namely the z axis, 
the rate equation for the n-phonon occupation probabil- 
ity p n is GEIEE] 



dt 



-Pn 



{ (n z + l)A_P nz+1 
(n z + 1) A + + n z A_ 



Pn 



+ n z A. 



-Pn.-l} 



(Bl) 



where A- and A + are the cooling and the heating rate, 
respectively. We use the rate equations for the good- 
cavity regime in |32j where the laser detuning A is much 
larger than the total dipole decay rate 7 = T/2, the ion- 
cavity coupling strength g and the cavity field decay rate 

K. 

In [35], the rates are derived presuming a stationary 
state solution. However, when the cooling is only along 
the axial direction, the radial phonon occupation number 
(n r ) increases over time due to recoil heating, 



dt 



O 2 

2 A 2 ' r 



(B2) 



Consequently, as seen in Eq. (A12|, the cooperativity 



decreases over time and results in a smaller photon scat- 
tering rate into the cavity. This eventually leads to a 
time dependence of the cooling and heating rates. Thus, 
A± rates never stabilize to stationary values due to their 
dependence on (n r ) and no steady-state temperature is 
achieved. 

Although the stationary state condition is not fulfilled 
here, we assume a pseudo steady-state which depends 
on the time-varying thermally averaged cooperativity. 
Thus, the resultant rates vary according to the thermally 
averaged cooperativity, 



Q 2 

± 7 2A 2 



<Pi 



C 



where a is a geometric factor (1/3 for Z' x , 2 — > J1/2 tran- 
sition), ipi is the cosine of the angle between the driving 
laser and z axis and S c is the relative detuning between 
the laser and the cavity resonance. 

In addition, we need to account for the cavity birefrin- 
gence. Because of our limited optical access, we can only 
probe at an angle of 45° to the vertical (y). Thus the 
probe couples equally to both birefringence modes of the 
cavity. In our case the cavity birefringence is not well 
resolved, and the best cooling rate is obtained when the 
cooling laser is tuned such that 5 C = — uj z ± <5&, where 
25b is the separation between the modes. In this case, 
the heating rate is not strongly affected by the birefrin- 
gence and is approximately the same as that in Eq. ( B3 1 , 
namely 



-4, 
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As for the cooling rate, the substitution of S c 
Eq. p3| leads to 
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into 



(B5) 



which translates into a lower cooling rate compared to a 
non-birefringent cavity. 



Multiplying Eq. (Bl I by n z and summing over all n z > 

(B6) 



results in a rate equation for (n 

d , , 



dt 



A_(n z )-A + ((n z ) + l) 



An analytical solution for (n z ) can be obtained by 
assuming a stationary state condition where A + P n _, — 
A_P„, + i for t — > 00 [32, 50 . Unfortunately, this condi- 
tion is not achievable with our current setup. Neverthe- 



less, even without an analytical solution, both Eq. (B2) 



and Eq. ( B6 ) are solved jointly and numerically and used 



(B3) for fitting to the experimental data. 
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